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ABSTRACT 



' Aims. We investigate the nature of the innermost regions with radii of several AUs of seven circumstellar disks around pre-main- 

I sequence stars, T Tauri stars in particular. Our object sample contains disks apparently at various stages of their evolution. Both single 

• , stars and spatially resolved binaries are considered. In particular, we search for inner disk gaps as proposed for several young stellar 
i-C ■ objects (YSOs). When analyzing the underlying dust population in the atmosphere of circumstellar disks, the shape of the lOfim 

feature should additionally be investigated. 

Methods. We performed interferometric observations in N band (8 - I3yum) with the Mid-Infrared Interferometric Instrument (MIDI) 
, at the Very Large Telescope Interferometer (VLTI) using baseline lengths of between 54 m and 127 m. The data analysis is based on 

* ■ radiative-transfer simulations using the Monte Carlo code MC3D by modeling simultaneously the spectral energy distribution (SED), 
^ ' N band spectra, and interferometric visibilities. Correlated and uncorrelated N band spectra are compared to investigate the radial 

distribution of the dust composition of the disk atmosphere. 

Results. Spatially resolved mid-infrared (MIR) emission was detected in all objects. For four objects (DRTau, RU Lup, S CrAN, and 
S CrA S), the observed N band visibilities and corresponding SEDs could be simultaneously simulated using a parameterized active 
^ ' disk-model. For the more evolved objects of our sample, HD 72106 and HBC639, a purely passive disk-model provides the closest 

1/^ I fit. The visibilities inferred for the source RULup allow the presence of an inner disk gap. For the YSO GWOri, one of two visibility 

, measurements could not be simulated by our modeling approach. All uncorrelated spectra reveal the 10/jm silicate emission feature, 

l/^ ■ In contrast to this, some correlated spectra of the observations of the more evolved objects do not show this feature, indicating a lack 

' of small silicates in the inner versus the outer regions of these disks. We conclude from this observational result that more evolved 

I dust grains can be found in the more central disk regions. 

' Key words. Infrared: stars - Accretion disks - Planetary systems: protoplanetary disks - Astrochemistry - Instrumentation: interfer- 

. ometer - Radiative transfer 

o 

> ■ 

•'H , 1. Introduction disks. Even the presence of a (stellar) companion and its location 

T . • . , , r, ■ , , ■ „ , , can be investigated using MIDI. The analysis of (interferomet- 

J, In this study we present the results of highly spatially resolved ^.^^ observations of more than one YSO at different evolutionary 
observations of seven pre-main-sequence stars acquired with the -^^ opportunity of verifying and constraining disk 

instrument MIDI at the VLTI. Our study is focused on YSOs, ^^^lution and planet-formation models (e.g., Leinert et al. 2004). 
TTauri stars in particular, because this class of object is consid- 

ered to represent the progenitors of solar-type stars. Studying the ^IDl, operating in the atmospheric N band, is sensitive to 
evolution of the disks around these stars with MIDI will provide ^^^^ ^>^qq K) matter, dust in particular, located in the inner 
us with an insight into the innermost and most dynamical regions of young circumstellar disks. It is generally assumed 

where the earliest stages of the formation of a planetary system ^^^^ ^^e formation of planetary systems starts in these inner disk 
can be found. regions (Klahr |2004t Wiinsch et aL |2005l Nagasawa et al. 2006J. 

In terms of spatial resolution, sensitivity, and spectral cov- ^ consequence of the planet-formation process, the disk struc- 
erage, MIDI is the only instrument suitable to the simultane- composition of the dust population is assumed to 

ous investigation of the temperature and density distributions in g^^j^g (Lisg^uer [19931 Gail 12003] Wolf i2007j). The dust com- 
the innermost disk regions, accretion effects, and dust proper- ^^^^^^^^ ^^^^ity distribution of old T Tauri objects with an 
ties. I.e., in determining th e evolutionary status of circumstellar > iq Myrs are assumed to have been strongly modified. In 

Send offprint requests to: A. A. Schegerer, e-mail: these systems, the accretion onto the£ar is sti'ongly reduced or 
schegererOmpia-hd. mpg.de 1^^^ stopped (Haitmann et al. |1993l ). Furthermore, the for- 

* Based on observations made with Telescopes of the European m^r circumstellar envelope of Cla ss I ob jects as well as the in- 
Organisation for Astronomical Research in the Southern Hemisphere nermost disk regions (Calvet et al. |2002| l may have disappeared 
(ESO) at the Paranal Observatory, Chile, under the programs 074.C- in these objects. In contrast, in YSOs with an age of ^1 Myr, the 
0342(A), 075.C-0064(A,B), 075.C-04I3(A,B), and 076.C-0356(A). accretion process usually plays a dominant role. 
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Binarity is a common phenomenon in YSOs. Extensive stud- 
ies in the near-infrared (NIR) wavelength range, where a spatial 
resolution of between 0.005" and 10" was reached, showed that 
between 23% and 55% of YSOs are binaries or even multiple 
systems (Simon et al. |1995t Kohler & Leinert 1998, Duchene et 
al. 1 19991 Ratzka et al. l2005l l. In a theoretical study, Artymowicz 
& Lubow d 19941 1 highlighted modifications to the structure of 
circumstellar and circumbinary disks resulting from a close com- 
panion. Furthermore, the nature of the infrared companions of 
TTauri stars remains a matter of debate in terms of the struc- 
ture and orientation of their circumstellar disks as well as their 
evolutionary stage (Koresko et al. 1997; McCabe et al. 2006). 

Simultaneous modeling of the spectral energy distribution 
(SED) - averaged over the entire disk - and of interferomet- 
ric measurements - sensitive to the inner, warm regions - is the 
main goal of this study. Modeling the SED alone does not allow 
us to constrain unambiguously the disk density and temperature 
structure (e.g., Thamm et al. 1 19941 1 . Spatially resolved observa- 
tions, such as interferometric observations with MIDI are needed 
to reduce the ambiguity caused by considering solely the SED. 
The modeling approach should be as simple as possible to re- 
duce the mathematical underdetermination of the model. A key 
issue of this study is to determine whether all objects in our het- 
erogeneous sample can be modeled with one single approach or 
if modifications and extensions of the model must also be con- 
sidered. Finally, it has not yet been shown whether any quantities 
of the objects, such as the stellar luminosity, stellar temperature, 
and mass accretion rate derived in previous, spatially unresolved 
observations can also be confirmed with interferometric, small- 
scale (range of AU) observations. 

Any deviations in the outcome of a homogeneous disk model 
from the (spectral and spatial) measurements point to (small- 
scale) perturbations of the disk structure. A non-representative 
selection of potential perturbations are an inner-disk gap (Calvet 
et al. 120021 Mathieu et al. il99Tl l, vorti ces (K lahri2004), a puffed- 
up inner-rim wall (Dullemond et al. I2001I I. and a close stellar 
companion (Mathieu et al. 119971 Artymowicz & Lubow [T994| l. 
Non-homologous evolution of the disk produces these perturba- 
tions (Wood et al. 2002, McCabe et al. l2006] l. 

MIDI provides both correlated, i.e., coherent spectra sen- 
sitive to the inner part, and uncorrected spectra sensitive 
to regions that can be non-interferometrically resolved. For 
the TTauri objects TWHya (Ratzka et al. |20U7] l and RYTau 
(Schegerer et al. 2008), an increase in the relative crystalline 
mass contribution to the lOfim silicate feature could be deter- 
mined for an increasing spatial resolution obtained with MIDI. 
This result was even found for the active object FUOri (Quanz 
et al. 1200 5). This finding was explained by the hypothesis that 
silicate dust evolves more rapidly in inner disk regions, i.e., at 
small radii r, than in outer regions of the disk (e.g., van Boekel 
et al. 1200 4*) assuming similar initial dust compositions. It is not 
known, so far, if this increase occurs only for intermediate-mass 
YSOs such as RYTau or if such a coiTelation exists even for 
fainter (~1 L©), younger (~1 Myrs old), and less active TTauri 
objects. Radial mixing could reverse the increase in evolved dust 
towards the more inner regions (Gail 120031) . In this study, we 
compare silicate features that effectively arise at different disk 
radii within an object. By focusing on one single object, the re- 
sults of this investigation do not depend on disk properties such 
as the evolutionary stage, initial composition, and disk inclina- 
tion but only on the origin of the feature. 

This paper is structured in the following way. In Sect. |2l we 
present the object sample and briefly describe our observations 
and data reduction. In Sect. [3] our modeling approach is outlined. 



The derived individual disk parameters are presented in Sect. |4] 
An analysis of the 10 //m silicate feature is given in Sect.|5] fol- 
lowed by a summary in Sect.|6] The Appendix Usts some further 
properties of the objects in our sample. 

2. Interferometric observations and data reduction 

Most of our objects are well-known TTauri stars. They belong 
to various star-forming regions such as the Taurus-Auriga (e.g., 
DRTau) and p Ophiuchi (e.g., HBC639) region. Selected ob- 
ject parameters are compiled in Table [T] (see Appendix |A] for 
more details). Apart from one, all objects are between 0.1 and 
3 million years old. The exception is HD 72106 B with an age of 
10 million years. Highly spatially resolved measurements in the 
NIR and mid-infrared (MIR) wavelength range have shown that 
the objects DR Tau and RU Lup do not have any stellar compan- 
ion (e.g., Ghez et al. [T997] l. In contrast, HD 72106 and SCrA 
at least are well-known binary systems. Considering the stellar 
mass, GWOri with - 3.7 M© would be the only Herbig 
Ae/Be object (M* > 3 Mq) in our sample, while in terms of 
spectral type only HD 72106 B would be a Herbig Ae/Be object 
(Vieira et al. |20U^ . 

The observing time with MIDI, the projected baseline length 
(B), and the corresponding position angle (PA) of the interfer- 
ometer during observations are listed in Table |2l A detailed de- 
scription of an observing sequence with MIDI is given in Leinert 
et al. (I2004I I. The calibration was completed using data for all 
calibrator stars that were observed during one night with the 
same observational mode as the scientific target and show clear 
visibility signals. The error bars in the visibility curves, which 
are presented below, represent the standard deviation when us- 
ing the ensemble of calibrator measurements. Since the sample 
contains various objects with known close-by infrared compan- 
ions, we checked whether these companions affect the visibility 
measurements. This examination was necessary because of the 
relatively wide field-of-view of MIDI, of about 1 ". The compan- 
ions do not affect the measurements. We used the MIA package 
(MIDI Interactive Analysis Software) for data reduction (Leinert 
et al. 2004; Ratzka 2005). The results were confirmed by using 
the independent EWS software. Both reduction software pack- 
ages are publicly available^ Several target/calibrator pairs were 
observed at different airmasses (Table |2]). Therefore, a coiTec- 
tion to the airmass was performed by multiplying the uncorre- 
cted spectra with the inverse of the cosine of the zenith angle. 
The visibility spectra do not need to be corrected because the 
correlated flux is normalized by the uncorrected flux. 

3. IVIodeling approach 

In the data analysis, we used the radiative-transfer code MC3D, 
which is based on a Monte Carlo method (Wolf et al. 119991 ) and 
was more recently developed to its present form (Schegerer et 
al. l2008l l. A predefined, parameterized model of a passive disk 
is the basic ingredient of our modeling approach. We use the disk 
model of Shakura & Sunyaev (.1973) , where the density distri- 
bution is given as: 



p{r,z) =pol — 1 exp 



z 



(1) 



available 



at 



' software packages are 

http : //www ■ itipia-hd . mpg . de/MIDISOFT/ 
and 

http : //www ■ strw . leidenuniv . nl/~ koehler/MIA+EWS-Manual/ 1 
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Table 1. Object properties derived in previous measurements or in this study. Different models of an object are provided in 
roman digits. Here, the coordinates (RA and DEC in J2000.0), distance (d), visual extinction (Ay), spectral type (SpTyp), 
stellar mass (M*), effective stellar temperature (T*), stellar luminosity (L*), and age are listed. Footnotes are references 
and point to the previous measurements. The footnotes "det." and "ass." represent parameter values that were determined 
and assumed in this study, respectively. 



Object (no.) 


RA [h m s] 


DEC [° "'] 


d [pc] 


Av [mag] 


SpTyp 


M, [Mo] 


T. [K] 


[Lo] 


Age [Myr] 


DRTau(I) 
















l-7(dct.) 




DRTau(II) 


04 47 06.2(2) 


+ 16 58 43(2) 


140(,,,, 


1.6(3) 


K7(3,4) 


0.80(4) 


4000(4) 


0.9(3) 


3(5) 


DRTau(III) 
















0-9(3) 




GWOri(I) 
GWOri (II) 


05 29 08.4(1) 


+ 11 52 13(1) 


440(6) 


1-3(6) 


G0(6) 


3.7(6) 


6000(6) 


40(de,.) 

62(6) 


1(6) 


HD 72106 B 


08 29 34.9(1) 


-38 36 21(1) 


290(5) 


0.0(7) 


A0(,) 


1-8(8) 


9500(7) 


28(det.) 


10(8) 


RULup 


15 56 42.3(1) 


-37 49 16(1) 


127(11) 


0.50(11) 


K8(ii) 


0-80(11) 


4000(11) 


1.3(11) 


1(11) 


HBC 639 


16 26 23.4(2) 


-24 20 60(2) 


170(11) 


5.7(11) 


KO(ii) 


2-0(11) 


4800(d„.) 


8.5(11) 


2(11) 


SCrAN(I,II) 
SCrAS 


19 01 08.6(2) 


-36 57 20(2) 


130(12,13) 


1-5(12) 
1.0(12) 


K3(i2) 
M0(i2) 


1-5(12) 

0-60(12) 


4400(12) 
3800(12) 


2.3(12) 

1.0(dct.) 


3(12) 
1(12) 


References: 1: Perrvman 09971; 2: 2 MASS catalogue (Cutri et al.l20031; 3: 


MuzeroUe et al. 


120031: 4: Mohanty et al. ( 2005 i; 5: Greaves 


120041: 6: 


Calvet et al. 120041: 



7: Schiitz et al. {20051 : 8: Wade et al. {20051 : 9: Vieira et al. {20031 : 10: Folsom (2007}; 11: Gras-Velazquez & Ray (2005) : 12: Prato et al. {20031 : 13: Johns-KruU et 
al. {20001 



The quantities r and z are the radius measured from the disk cen- 
ter and the vertical distance from the disk midplane, respectively. 
The quantity is the stellar radius. The function h(r) represents 
the scale-heightQ 

h(r) = hioo ( — f (2) 

with hiQo - h{r - 100 AU). The quantity po allows us to scale 
the total disk mass Mjisk. We assume a gas-to-dust mass ratio of 
100: 1 . In our approach the exponents a and j3 satisfy the follow- 
ing relation: 

« = 3(;6-i). (3) 

The latter relation results from the coupling of surface density 
and temperature of the disk (Shakura & Sunyaev[1973]l. We note 
that the approach established by Shakura & Sunyaev is correct 
only for geometrically thin disks, i.e., where z <K r. However, 
since the density exponentially decreases with increasing z and 
the scale-height h{r) is at least one order of magnitude smaller 
than the corresponding radius r, this approximation should be 
justified for the predominant mass fraction. The error in the den- 
sity of regions z ^ r is of the order z^/r^ for Taylor's series. The 
approach of Shakura & Sunyaev was already successfully used 
in modeling HH 30 IRS (Burrows et al. T996^, Wood e tal. [19981 
Cotera et al. 2001 ) and the Butterfly Star (Wolf et al. l2003l) . In 
contrast to the more basic, but time-consuming, approach of hy- 
drostatical equilibrium (Schegerer et al. 2 008) , the predefined, 
parametrized disk model of Shakura & Sunyaev ( 119731 ) avoids 
the difficulties in handling regions of high optical depth (e.g., 
Sonnhalter et al. 1993 ). However, even in predefined density dis- 
tributions, the resulting temperature distribution must converge 
independently of the number of photon packages emitted. 

For all objects besides HD72106B and HBC 639, accre- 
tion effects were also considered. The potential energy of a 
particle on its way towards the star is partly released in the 

^ The scale-height is defined as the vertical distance from the mid- 
plane, where the density has decreased by a factor e x 2.718 (Euler's 
constant). 



disk midplane assuming the canonical modeling approach of 
an geometrically thin, active disk as formerly established by 
Lynden-Bell & Pringle (|1974t , and Pringle ( |19811 l. Considering 
the "magnetically mediated" modeling approach (e.g., Uchida 
& Shibata [19841 Bertout et al. [19881 Calvet & Gullbring[ T998l l, 
the accretion disk is truncated by the stellar magnetic field at a 
radius /?bnd- Thereupon, more than half of the potential energy 
of the accreting particles is released in a boundary region above 
the stellar surface. This boundary region is heated to a temper- 
ature of Tbnd by the accreting material. The effects of the this 
approach were described and applied by Schegerer et al. (2008)1 
in modeling the circumstellar disk around RYTau. Calvet & 
Gullbring (1 19981 1 theoretically derived reasonable ranges for the 
latter both accretion values, i.e., 5700 K < Tbnd < 8800 K and 
2Ri, < 7?bnd < 5Ri,. We note that the effects of any variations 
in the boundary temperature Tbnd and the magnetic truncation 
radius /?bnd on the SED and MIR visibilities can be neglected in 
the above-mentioned ranges (D' Alessio et al. |1998l l. In our mod- 
eling approach, we always used standard values for the bound- 
ary temperature and truncation radius, i.e., Tbnd = 8000 K and 
Rbnd = 5 R*. These values were applied by Akeson et al. (2005]l 
and Schegerer et al. ( I2008I I for other T Tauri stars. The accretion 
rate M, which determines the total accretion luminosity Lace, is 
not an independent model parameter. But its value is constrained 
by the results of previous measurements. Correspondingly, the 
stellar mass M*, the effective stellar temperature T*, and the 
stellar luminosity L*, which were derived in previous studies 
were starting parameters in the modeling. In some cases, it was 
necessary to deviate from these results to be able to simulate 
the measurements considered in this study. These deviations are 
discussed below. 

In our modeling approach, we assumed a dust mixture of 
"astronomical silicate" and graphite with relative abundances of 
62.5% and 37.5%, respectively (Draine & MalhotralT993). We 
used a two-layer disk model, i.e., the disk interior, where the 
optical depth in N band fulfills the condition tn > 1, contain- 
ing dust grains with a maximum particle size of fln,ax - 1 mm. 
A maximum particle size of Omax = 1 mm was found for sev- 
eral T Tauri stars with millimeter observations using the Very 
Large Array (Rodmann et al. |2006l l. Here, the optical depth tn is 
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Table 2. Overview of our observations with MIDI. The date and the observing time (UT) as well as telescope pairs {TP), 
and both the projected (effective) length (B) and angle of the baseline {PA, measured from North to East) during the 
observations are presented for targets and corresponding calibrator stars. The given airmass (AM) is an average value for 
the observation. In the following columns, we add the diameter of the calibrator star and the coiTesponding references. 
The diameters were derived from fitting the SED of the objects. In the last column, we also provide comments concerning 
the observation and data reduction. 



Date 




f IT 




vjDject 


I r 


B [mj 


tA L J 


A.M 


Diameter [mas] 


Kei. 


Comments 


Jan. Ul zUUj 


z : 


jU — 2 : 


59 


DRTau 


U 13-U 14 


C 1 

oi 


1 AA 

luo 


1 1 






*'**** 


Jan. Ul ZUUj 


Q . 

D . 


Ij — J . 


1 n 
1 / 


1 I FA Q 1 /I O 1 


U 1 J-U 14 


An 


1 riA 


I. J 


Z. /4 ± U.Uj 


1 
1 




Jan. Jl zUUj 


n . 

J . 


JD — J . 


JO 


\JT\ '2 1 /I O 1 


U 1 J-U 14 


JO 


1 n/1 
1U4 


1 'J 
I. J 


z. /4 ± U.Uj 


1 
1 




Mar OP^ 9nnS 






• j4 


1 1 L» J i '-rZ. i 


UT3-UT4 




105 






\ 




Mar. 01" 2005 


: 


21 - : 


35 


GWOri 


UT3-UT4 


56 


105 


1.3 








Mar 03"^ 9005 


: 


02 - : 


13 


HD 3 1 47 1 


UT2-UT4 


87 


79 


1.3 




\ 




Mar. 03"' 2005 


: 


28 - : 


45 


GWOri 


UT2-UT4 


88 


79 


1.3 






* 


Dec. 3r' 2005 


4 : 


19-4 : 


39 


HD 72106 


UT1-UT4 


127 


41 


1.1 






** 


Dec. 3P' 2005 


4 : 


53-4: 


56 


HD 69142 


UT1-UT4 


130 


48 


1.1 


2.18 ±0.18 


2 




Mar. 11"' 2005 


: 


29 - 1 : 


01 


HD 72106 


UT2-UT4 


88 


71 


1.1 






** 


Mar. 11"" 2005 


1 : 


08 - 1 : 


27 


HD 69142 


UT2-UT4 


89 


78 


1.0 


2.18 ±0.18 


2 




Aug. 26^^' 2005 


: 


52 - : 


56 


HD 152885 


UT2-UT4 


84 


96 


1.1 


2.88 ± 0.09 


2 




Aug. le' 2005 


1 : 


55 - 2 : 


25 


RULup 


UT2-UT4 


64 


122 


1.5 






* 


Aug. 2005 


2 : 


38-2 : 


40 


HD 152885 


UT2-UT4 


69 


116 


1.4 


2.88 ± 0.09 


2 




Aug. 26" 2005 


3 : 


02 - 3 : 


04 


HD 178345 


UT2-UT4 


84 


96 


1.1 


2.42 ± 0.03 


2 




Aug. lb /UUj 


A ■ 


AQ A • 

46 — 4 . 


jU 


rlU 1 /?5J4j 


T TTO T ^TA 

U Iz-U 14 




1 1 A 

1 10 


1 A 

1.4 


LAL ± U.Uj 


L 




May 15*2006 


2 : 


29 - 2 : 


43 


RULup 


UT2-UT3 


46 


17 


1.2 








May 15*2006 


4 : 


51 -5 : 


03 


HD 152820 


UT2-UT3 


46 


13 


1.3 


2.63 ±0.15 


2 




May 25* 2005 


3 : 


20 - 3 : 


28 


HD 152820 


UT3-UT4 


55 


94 


1.1 


2.63 ±0.15 


2 




May 25" 2005 


3 : 


43 - 3 : 


50 


RULup 


UT3-UT4 


61 


103 


1.0 








May 25" 2005 


7 : 


32 - 7 : 


40 


HD 152820 


UT3-UT4 


58 


130 


1.1 


2.63 ±0.15 


2 




May 25"' 2005 


7 : 


55 - 8 : 


02 


RULup 


UT3-UT4 


55 


149 


1.5 








Apr. 18'" 2005 


4 : 


11 -4 : 


18 


HD 142804 


UT2-UT4 


71 


66 


1.3 


2.80 ± 0.08 


2 




Apr. 18"' 2005 


4 : 


39 - 4 : 


58 


HBC 635 


UT2-UT4 


75 


61 


1.2 








Anr 19"'2005 


3 : 


19 - 4 : 


33 


HD 142804 


UT2-UT4 


60 


58 


1.5 


2 80 + 08 


2 




Apr. 19"' 2005 


3 : 


44 - 4 : 


11 


HBC 635 


UT2-UT4 


66 


52 


1.4 






* 


Aug. 25"" 2005 


23 : 


19-23 


: 21 


HD 142804 


UT2-UT4 


60 


116 


1.0 


2.80 ± 0.08 


2 




Aug. 25'" 2005 


23 : 


45-23 


: 46 


HD 142804 


UT2-UT4 


58 


118 


1.0 


2.80 ± 0.08 


2 




Aug. 26'" 2005 


: 


08-0: 


27 


HBC 635 


UT2-UT4 


59 


121 


1.1 








May 31" 2005 


10 : 


07 - 10 


: 38 


SCrAN 


UT2-UT4 


70 


109 


1.4 






=j= * * 


May 31" 2005 


10 : 


48 - 10 


: 49 


HD 178345 


UT2-UT4 


66 


120 


1.5 


2.42 ± 0.03 


1 




Jun. 28'" 2005 


7 : 


57-7 : 


58 


HD 178345 


UT3-UT4 


59 


137 


1.2 


2.42 ± 0.03 


1 




Jun. 28'" 2005 


8 : 


17 - 8 : 


39 


SCrAN 


UT3-UT4 


56 


145 


1.4 








Aug. 26'" 2005 


: 


52-0 : 


57 


HD 152885 


UT2-UT4 


84 


96 


1.1 


2.88 ± 0.09 


2 




Aug. 26'" 2005 


2 : 


38-2 : 


40 


HD 152885 


UT2-UT4 


69 


116 


1.4 


2.88 ± 0.09 


2 




Aug. 26'" 2005 


3 : 


02 - 3 : 


04 


HD 178345 


UT2-UT4 


84 


96 


1.1 


2.42 ± 0.03 


1 




Aug. 26'" 2005 


4 : 


04 - 4 : 


35 


SCrAS 


UT2-UT4 


72 


110 


1.3 






* 


Aug. 26'" 2005 


4 : 


48 - 4 : 


50 


HD 178345 


UT2-UT4 


69 


116 


1.4 


2.42 ± 0.03 


1 





* bad weather conditions (not perfect beam overlap, noisy data) 
** data reduced using EWS package, only 

* * * results from MIA and EWS package are not consistent at short wavelength (deviation: 13%) 

* * * * results from MIA and EWS package are not consistent for the whole band (devation: 14%) 
References: I : Cohen et al. U999) : 2: MIDI consortium (in prep.) 



measured for constant radii, from the disk atmosphere vertically 
to the disk midplane. The surface regions consist of interstellar, 
non-evolved dust with amax — 0.25 yccs\ as found in the interstel- 
lar medium (Mathis et al.'1977i In both the disk interior and the 
disk surface layer a minimum grain size of Omin - 0.005 /im and 
a grain size distribution of cC^ '^ for a e [amin, ^max] was adopted 
(MRN grain-size distribution; Mathis et al. 1 19771 1. 

After the temperature distribution was determined assuming 
a density given by Eq.[T] the SED, the projected image of the star, 
and its circumstellar environment, assuming an inclination angle 
were calculated. In the modeled images from which the visi- 
bilities result, we generally used a spatial resolution of <0.2 AU 



in the MIR range, which is a factor ~10 higher than the spa- 
tial resolution reached in our MIDI observations. For the images 
from which we derived the NIR visibilities, we used a resolu- 
tion of <0.03 AU, i.e., a factor of 10 higher than the resolution 
reached by any existing NIR interferometer. 

Independent model parameters of our approach are the outer 
radius 7?out, the scale-height /iioo, the exponent/?, and the inclina- 
tion angle of the circumstellar diskjj Assuming a certain dust 
mixture in the models, the disk mass is fixed predominantly by 
the flux in the millimeter wavelength range. Considering a stel- 



^ An inclination angle of 0° corresponds to a face-on disk. 
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lar temperature T*, the quadratic distance law, a mean dust sub- 
limation temperature of 1500 K, and the specific absorption co- 
efficients K of the adopted dust set, the temperature of single dust 
grains at each radius can be determined. However, our approach 
is only valid for optically thin media. For optically thick me- 
dia, the grain temperature at the inner-disk edge is even higher 
because back radiation from outer dust grains produces an ad- 
ditional heating of the grains at the inner edge. The resulting 
sublimation radius R^uh is the initial value of Ri^ in our modeling 
approach. 

Initially, the simulation grid of independent parameters is 
scanned coarsely during modeling. Simultaneously, we consider 
parameter values that were derived in preceding studies. By 
means of manual modifications of the model parameters, we 
search for the model that reproduces most successfully the mea- 
sured SED and MIR visibilities, simultaneously. The SED and 
the MIR visibilities that we try to model are complex functions 
of many diflFerent disk and dust parameters. However, the effects 
of any modification of the parameters on modeled SED and vis- 
ibilities can be estimated. A manual modification is necessary 
because a disk model with a specific SED and images at differ- 
ent wavelength^ takes about an hour to compute. In our search 
for the best-fit model the uniqueness of our final models cannot 
therefore be proven. However, we verified whether the modeling 
results can be improved by varying the model parameters using 
the following step widths: 

AL* = 1 Lq, AMdisk = 0.5 Mdisk, 
A7?out = 10 AU, A/Jin = 0.05 AU, 
A;8 = 0.1,and A/iioo = lAU. 

Considering the error bars of the measurements, the modeling 
results could not be improved using finer step widths. The step 
widths can be considered as the precision to which the local min- 
imum in the;t'^-surface can be determined. The determination of 
the modeling errors would require a fit of the ;t'^-surface in the 
simulation grid to the independent parameters around the deter- 
mined local minimum. 

4. Results 

All scientific targets could be well resolved, i.e., the visibili- 
ties are lower than unity. Table [3] presents the parameters re- 
sulting from our effort to model the circumstellar disks of seven 
YSOs. In the following, we present and discuss our modeling 
results for each object in our sample, in detail. Previous obser- 
vations, including photometric measurements, are mentioned in 
AppendixlAl 

4.1. DRTau 

Considering the low gradient of the SED at wavelengths A ^ 
1 yum (Fig.lTJ, it is evident that the (measured) SED of DRTau 
is represented by a black body with the stellar temperature of 
r* - 4050 K and an additional radiation source. The accre- 
tion process is the other strong influence on the radiation emit- 
ted in this wavelength range. In addition to the stellar luminos- 
ity and the accretion rate, the parameters of the model (I) of 
our approach confirm the values derived in the previous study 
of Akeson et al. (2005 ) using long-baseline NIR interferometry 
from the Palomar Testbed Interferometer (PTI). However, they 
found a stellar luminosity of = 0.9 Lq, which is a factor of 
~2 lower than the value found in our model. Simultaneously, 

* used for the calculation of the modeled visibility 



the accretion rate of their model, i.e., M = 8 x lO^^Moyr"', is 
four times higher than the value determined in our first model. 
However, SED and MIR visibilities can also be reproduced by 
a model (II) by considering the stellar luminosity and the ac- 
cretion rate found by Akeson et al. (2005) while keeping all 
other model parameters constant. Both models similarly fit MIR 
visibilities and SED to within the error bars. We conclude that 
the intrinsic stellar luminosity and the accretion luminosity can- 
not be disentangled considering SED and MIR visibilities, only. 
The derived accretion luminosity Lace depends on the applied 
accretion rate M: the model (I) has an accretion luminosity of 
Lace = 0.19 Lq for M - 2x 10'^ Mgyr"' while an accretion lu- 
minosity of Lace = 1 .0 Lq for M = 8 X 10"** Mgyr"' is emitted by 
the model (II). The sum of the intrinsic stellar luminosity L*^ and 
the accretion luminosity Lace in both models equals ~1.9Lsun- 
Finally, we mention that Mohanty et al. (2005) derived an ac- 
cretion rate of M = 2.2 x lO^^Moyr"' analysing Call lines. 
Robitaille et al. (120071) determined a lower and upper limit to 
M, i.e., 6.9 X lO-'^Moyr"' ^ M ^ 1.0 x lO^'^Moyr"' assum- 
ing the SED, only. If the NIR visibility data obtained by Akeson 
et al. (I2005I I is simultaneously taken into account, we find that 
the disks of models (I) and (II) appear too spatially resolved, re- 
sulting in too low NIR visibilities. Only an additional reduction 
in the inner disk radius to = 0.05 AU improves our fit to 
the NIR visibility data in the model (II). A similar improvement 
cannot be reached in model (I) by reducing its inner disk radius. 
Our model of DRTau with L* = 0.9 Lq, M = 8 x 10 ** Mgyr 
and Ri„ - 0.05 AU is called model (III). The lower right panel in 
Fig. [1] shows measurements and model output for the NIR vis- 
ibility of model (III). We note that the fit to the MIR visibility 
data at wavelengths of 12.5 /im becomes poorer in model (III) 
and deviates from the measured data by < 7% as the gradient 
of the modeled MIR visibility curve decreases. We summarize 
that model (III) is the only model that reproduces all available 
data sets apart from slight deviations in the MIR visibilities at 
~12.5 fim. 

As we mentioned in Sect.|3] the sublimation radius is approx- 
imated considering a stellar temperature of L* = 4000 K and the 
absorption coefficients of the dust composition that we used in 
our two-layer disk model. This approach is only an approxima- 
tion because we assume an optically thin disk in determining 
of the sublimation radius /?;„. The true sublimation radius Ri„ 
should be larger because the back radiation of adjacent grains 
is neglected. We obtain a sublimation radius of 7?sub ~ 0.06 AU 
even for the dust for which Omax - 1 mm because of a low num- 
ber of large grains with a > 0.25 /im. Since the inner radius Rin of 
the modeled disk is slightly smaller than the sublimation radius, 
this value of ^sub could imply that larger grains exist in the disk, 
in the innermost disk in particular, than the MRN grain-size dis- 
tribution predicts (Mathis et al. 1977). This result is also found 
in Sect. |5] and predicted by Isella & Natta (2005). However, the 
hint for larger grains is only weak since the large error bars of 
the NIR data still allow a slightly larger inner-disk radius Ri^ 
than we have assumed in model (III). 

As the disk crosses the line of sight of the observer with in- 
creasing inclination angle, the NIR flux at ~2.2 //m of the model 
of DR Tau decreases by 38% in steps of A-d- - 5° for small incli- 
nations. Therefore, the inclination angle can be determined with 
high accuracy. If the measured visual extinction Ay ~ 1 .6 mag 
(Muzerolle et al. |2003j can be ascribed exclusively to the cir- 
cumstellar disk around DR Tau that crosses the line of sight of 
the observer, an inclination of 20° can be derived. If the ex- 
tinction is also caused by interstellar media outside the system. 
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Table 3. Parameters of the models that simulate the measured SED and MIR visibiUties best. Footnotes of 
the mass accretion-rate point to measurements in previous studies. Values that are marked with (det.) are 
determined in this study. 



object (model-no.) 


Mr t [Mol 


R , TAUl 






/7inn [AUl 
"l()u L-'^'^J 


L J 




DRTau(l) 






0.10 








2 • 10"* 


DRTau (II) 


0.1 


90 


0.10 


0.75 


15 


20 


8.0 • 10-* 


DRTau(III) 






0.05 








8-0 -10^1, 


GW Ori (I,II) 


1.0 
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0.35 


1.10 


22 


10 
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HD 72106 B 


0.005 


40 
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1.30 
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60 
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RULup 


0.1 
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0.10 


0.90 


20 


28 


1.0- lo-j* 


HBC 639 


0.1 
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0.10 


1.00 


10 


65 


NoAccretion 


SCrAN(I) 


0.03 
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0.05 


1.10 


9 


10 


4 • 10-* 


SCrAN(II) 








0.90 


12 


45 
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References: 1: Akeson et al. 120051 ; 2: Calvet et al. 120021 : 3: Herczeg et al. 120051 



which is not considered in the model, the derived inclination an- 
gle is correspondingly smaller 

Our observation of DR Tau with MIDI achieved a theoret- 
ical spatial resolution of % 17 mas (~2AU). To this limit, the 
visibility data do not provide any evidence of a companion 
with a brightness ratio close to unity, which would have pro- 
duced a characteristic sinusoidal signature in the visibility (e.g., 
Schegerer et al. |2008l l. We note that the detection of a close bi- 
nary is difficult when the binary is oriented nearly orthogonal to 
the baseline vector of our interferometric observations. 

Our finding confirms the results of former speckle- 
interferometric measurements in the NIR range (Ghez et 
al. 1 19971 1, where a theoretical spatial resolution of 53 mas could 
be achieved. Two previously unpublished lunar occultation data 
sets of high quality (SNR » 20 - 40) are available for DRTau 
from the 3.5 m telescope in Calar Alto, the first one recorded 
in November 1997 with an InSb fast photometer and the sec- 
ond one in September 1999 with the IR array OmegaCass in fast 
readout for a 32 pixel window. A detailed analysis allows us to 
conclude that the source appeared point-like with an upper limit 
of (2.70 + 0.65) mas, i.e., (0.4 + 0.1) AU and (5.0 + 0.4) mas, 
i.e., (0.7 + 0.06) AU in both data sets. Upper limits can also be 
set to the flux because of a possible circumstellar emission, and 
they are ^ 5% over 160 mas and ^2.5% over 420 mas, for both 
observations. 

4.2. GWOh 

As former millimeter measurements have shown (Mathieu et 
al. 119951 ) and can be confirmed by this study, the circumstel- 
lar disk around GWOri is very massive, i.e., Mdisk = 1.0 M©. 
Furthermore, compared to other objects in this study, the disk of 
GW Ori has the largest outer radius Rqm- Besides the luminosity, 
all other stellar parameters and the accretion rate found in this 
study are consistent with the result of Calvet et al. (2004), who 
analyzed the UV spectrum of this source. The accretion lumi- 
nosity is Lace = 5 Lq in our model. If an intrinsic stellar lumi- 
nosity of L* = 62 Lq is used, as was determined by Calvet et 
al. (I2004I I. the resulting flux in K band would exceed the pho- 
tometric measurement by 20%. In Fig. |2l the SED and MIR 
visibilities of a model with a stellar luminosity of L* = 62 L© 
is also shown (model (II): gray curves/visibility bars). Although 
both the SED and the visibilities for Z? = 56 m can be repro- 



duced, no model could be found that reproduces simultaneously 
the visibilities of a baseline of B = 88 m. For spectroscopic mea- 
surements, Mathieu et al. d 19971 1 found that GW Ori has a stellar 
companion at a projected separation flsep between cisep = 1.08 AU 
and flsep - 1.18 AU and of stellar mass between = 0.5 M© 
and Mi, - 1 Mq. According to Mathieu et al. ( 1997), this com- 
panion creates a dust-free ring in the circumstellar disk between 
r = 0.17 asep and r - 3.0asep. By means of a theoretical in- 
vestigation, Artymowicz & Lubow (1994) determined an inner 
radius of ~2.3 asep for the circumbinary disk of the system. They 
claimed that each component has its own circumstellar disk of 
outer radii of 0.46 asep and 0.20 asep. 

In our approach, the inner disk radius equals the sublima- 
tion radius {R\a - Raab)- We do not consider a dust- free in- 
ner gap as proposed by Mathieu et al. (|J_997]) and Artymowicz 
& Lubow d 1 9941 ). However, in another model a dust-free, ring- 
shaped gap is cut in our disk model. We therefore assumed the 
same disk parameters and /iioo for the circumbinary and cir- 
cumstellar disk. The disk gap extends from the radius 0.5 AU to 
3.0 AU. In the latter model, we obtained the following results 
(not shown in the figure): the NIR flux obtained from this model 
decreases by 70% and cannot fit the NIR range of the SED any- 
more. Although the modeled visibilities V{A) for the long base- 
line B = 87 m corresponds to a positive gradient, the visibilities 
for the short baseline B = 56 m decrease by 20% across the com- 
plete N band. Any modifications of the radii that limit the disk 
gap could not even improve our model considering the measure- 
ments. 

We assume that the disk gap is probably neither dust-free 
nor ring-shaped. In a theoretical study, it was shown that ma- 
terial streams between the circumbinary and the circumstellar 
disks guarantee the accretion onto each component (Giinther & 
KIey2Q02 ). Even the stellar companion that is not implemented 
in our approach could affect the MIR visibilities and SED. A 
strong concentration of MIR intensity in the innermost disk re- 
gions that could not be spatially resolved with MIDI and an 
abrupt truncation of the MIR intensity distribution potentially 
caused by a disk gap result in the increase of the MIR visibility 
on long baselines. 

To reveal the innermost disk region of GW Ori, more compli- 
cated modeling approaches as well as additional interferometric 
observations are essential. Measurements with the NIR interfer- 
ometer AMBER at the VLTI, for instance, could provide infor- 
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Fig. 1. SED and MIR visibilities for a projected baseline of B = 60 m obtained from the measurements and our models (I), (II), and 
(III) of DR Tau. The model parameters are listed in Table[T]and Table[3] Upper row: SED and an enlargement showing photometric 
data from the visible up to the NIR wavelength range. The curves are ascribed to the specific models by numbers. Real photometric 
data are plotted with eiTor bars. The dashed lines represent the intrinsic stellar flux of models (I), (II), and (III). Lower-left figure: 
The modeled visibilities were derived from the corresponding model images for the wavelengths of 8.5 jum, 9.5 yum, 10.6 jum, 
1 1.5 fim, and 12.5 yum. The colored bars represent intervals that limit the MIR visibilities V{A) for different position angles but the 
same inclination of the models (I: dark grey) and (III: light gray). For clarity, the visibilities obtained from model (II) are not shown 
but they are located between the results of model (I) and model (III). The measured data are included with error bars (black, solid 
line). Lower right figure: NIR visibilities as a function of the projected (effective) baseline length B of an interferometer derived for 
a wavelength of 2.2 fim that result from model (III: light gray). Model (III) is the only model that fit all available data sets. The dark 
grey line represents the output of model (I) and (II), respectively. The measured NIR visibihties obtained with the PTI (Akeson et 
al. l2005l l are included with error bars. 



mation about the small-scale structures of the innermost regions 
that potentially deviate from axial symmetry. 

4.3. HD 72106 B 

HD72106B is the least well-studied object in our sample. In 
our modeling approach, we use a luminosity of = 28 Lq. 
This value results directly from a fit by a Planck function of 
temperature - 9500 K (Schutz et al. 120051 ) to the photomet- 
ric measurements in U-, V- and R-band. Two properties make 
HD 72106B a special object in our sample: 

First, with an age of lOMyr, HD72106B is a relatively old 
TTauri object with weak infrared excess. The equivalent width 
of the Ha line that is smaller than 10 A (Vieira et al. |2003| l is a 
sign of weak accretion. Our model of HD 72106 B consists of a 
passive disk without considering any accretion effects. For sev- 



eral TTauri objects of similar age (e.g., TWHya; s. Calvet et 
al. '2002) an inner, dust-free gap in the range of ^ 1 AU, i.e., 
larger than the sublimation radius could be detected. However, 
because the SED and MIR visibilities of HD 72106 B can be sim- 
ulated by assuming that the inner edge of the (dust) disk equals 
the sublimation radius, a shift in the inner-disk edge to larger 
radii is not obvious. 



Second, HD72106B is the infrared companion of 
HD 72106 A with a projected distance of 280 AU. HD72106A 
has already developed into a main-sequence star (Wade et 
al. I2005I I. The main component is probably responsible for the 
truncation of the outer disk around the B component considering 
an outer radius of 40 AU. The small outer radius of the disk 
finally produces a strong decrease in the FIR flux (s. Fig|3]l. 
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Fig. 2. SED and MIR visibilities for the projected baselines of B - 56 m and B = 88 m obtained from the measurements and 
our model of GWOri. The results of two different models are plotted that differ in the stellar luminosity (dark gray: = 40 Lq, 
light gray: - 62 Lq). The increasing visibility of the longest baseline B - 87 m, which cannot be reproduced by our modeling 
approach, could point to a bright, truncated density structure in the inner disk region (< 5 AU), although different scenarios are also 
possible. 



4.4. RULup 

As Lamzin et al. (1996) has shown, RULup was at least tem- 
porarily a strongly accreting object. By analyzing its SED, they 
derived an accretion rate of M = 3 x 10"^ Moyr"'. In this con- 
text, Reipurth et al. (119961 ) found that the Ha emission has broad 
wings of up to 900 kms"' in width. Furthermore, the Ha emis- 
sion of RULup shows one of the broadest equivalent widths 
among TTauri stars (Giovannelli et al. ll995l l. However, by con- 
sidering SED and MIR visibilities, we can exclude such a high 
accretion rate M. We determine a maximum value of M = 1 x 
10"^ Moyr"' with an accretion luminosity of only Lace = 0.1 L©. 
Herczeg et al. (|2005i l found that M = (5 ± 2) x lO^^Moyi--' 
which supports our result. A higher accretion rate produces a 
decrease in the MIR visibility because the irradiation of the disk 
increases. Simultaneously, the NIR and MIR flux would exceed 
the photometric measurements. Figure |4] represents our best-fit 
model with M = 1 x lO^^Moyr"' (black curve/visibility bars). 
The stellar luminosity of = 1 .27 L© corresponds to a value 
derived by Gras- Velazquez & Ray (120051) . It represents the av- 
erage of previous measurements (Herczeg et al. '2005 : 0.49 L©; 
Niirnberger et al. 1997: 2.2 L©). A model that simulates the flux 
in the N band, i.e., the silicate feature, could not be found. The 
relative deviation in the MIR range is 37%. The visibilities ob- 



tained on the baselines B = 61 m and B = 65 m and for similar 
position angles significantly differ, which could be due to a mea- 
surement error. 

Based on spectro-astrometric observations, Takami et 
al. (120011 ) suggested that an inner gap in the circumstellar disk 
of RU Lup exists between ~2 AU and ~4 AU (20-30 mas) as jus- 
tified by the following arguments. Because of the positional dis- 
placement of the extended wings, the emission of forbidden lines 
such as [01] and [SII] in the spectra of RU Lup is assumed to 
arise in a collimated large-scale outflow. The apparent lack of 
redshifted line components results from an obscuration of the 
dorsal outflow by the circumstellar disk (Eisloffel et al. 120001 ). 
Takami et al. (2001) found that the Ha line in the spectrum of 
RU Lup, which is evident predominatly in the innermost regions, 
shows both blueshifted and redshifted wings. This finding is ex- 
plained by an inner disk gap that allows regions behind the disk 
to be observed. The gap could be evoked by an unseen (plan- 
etary) companion. Since the wings of the Ha line shows the 
same positional displacement as the forbidden lines, Takami et 
al. (120011) suggested that the Ha line originates predominantly 
in the outflow. 

Similar to our modeling approach of GWOri, we defined 
a dust-free gap in the disk model that is presented in Fig. |4] 
(Tables [TT3T ). The gap size was varied using values of Ri^ of up to 
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Fig. 3. SED and MIR visibilities for the projected baselines of B = 88 m and B - 111 m that results from the measurements and 
our model of HD72106B. Dotted curves represent remnants of the telluric ozone band that could not be eliminated by the data 
reduction. 



3 AU for the inner gap radius, and 3 AU and 5 AU for the outer 
gap radius. An outflow was not considered in our model. The 
resulting changes can be summarized as follows: as shown in 
Fig. El a gap lowers the NIR flux by -7%. We note that t he NIR 
flux is highly variable during a year (Giovannelli et al. I1995I ). 
Previous photometric measurements (Giovannelli et al. 1 19951 1 
would be consistent with lower NIR fluxes. A disk gap as used 
here would even lower the MIR flux and slightly improve our fit 
to the N band. The MIR visibility simultaneously decreases at 
wavelengths < 10.5 ^va as the relative MIR flux contribution of 
outer disk regions in relation increases. For longer wavelengths 
> 10.5 yum, the disk becomes more compact as the visibilities 
slightly increase. Our visibility measurements are inconsistent 
with a disk gap and an inner radius of r = only, but a gap 
of between 2 AU and 4 AU fits the visibility data as well as the 
model without gap. 



Although there are observational hints for a disk gap, we are 
unable to conclude whether the model with a disk gap of inner 
radius ^ 2 AU or a model disk without gap provide a closer fit to 
the data. Further interferometric observations and simultaneous 
photometric measurements in the NIR wavelength range would 
allow a final confirmation or rejection of the presence of a gap 
in the disk of RU Lup. 



4.5. HBC639 

Apart from HD 72106B, the source HBC639 is the only object 
where accretion has not been considered. Although HBC639 
has an age of between 1 and 3 million years (Gras-Velazquez 
& Ra y I2005l l. we can neglect the derived accretion rate (Prato 
et al. I2003I I. for which we estimate an upper limit of M < 
1 X 10~^Moyr Since the equivalence width of the Ha line 
is smaller than 10 A, HBC 639 is not a classical T Tauri star but 
belongs to the class of weak-line T Tauri stars. 

The parameter combination listed in Table[3]could reproduce 
the MIR visibilities. However, this model produces a strong de- 
cline in the FIR flux and is therefore unable to reproduce the pho- 
tometric upper limits in the FIR range. HBC 639 has an infrared 
companion at a projected distance of 320 AU (s. Appendix [All. 
The strong decline in the FIR flux could be explained by the 
presence of an infrared companion that truncates the outer-disk 
region of the main component. In the latter case, the measured 
flux in the FIR (s. Fig.|5]l can mainly be ascribed to this infrared 
companion. The companion is again a binary, as well, and deeply 
embedded in a circumbinary envelope. 

A larger value for the disk parameter yS results in a more 
flared disk. Outer disk regions can be more eff'ectively heated 
and the intensity distribution in the MIR range therefore de- 
creases less strongly. A consequence is a decrease in the MIR 
visibilities. Fig.|5]also shows the results of a model with p - l.l 
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Fig. 4. SED and MIR visibilities for the projected baselines of Z? = 46 m, Z? = 55 m, Z? = 61 m, and Z? = 64 m obtained from the 
measurements and our model of RULup. For comparison, the light gray and dark gray colored curves/bars are results from two 
different models with and without disk gap, respectively. The black region of the visibility bars represents an intersection. The gap 
-if present- stretches from 2 AU to 4 AU according to Takami et al. ( 1200 11 1. 



for which none of the other parameters were modified. We note 
that the visibilities decrease by at least 10% for all baselines and 
the FIR flux increases by almost a factor of 2. 



4.6. SCrAN 

Only the photometric measurements in the NIR and MIR range 
can exclusively be ascribed to the Northern component of this 



binary. All other photometric measurements could not spatially 
resolve the individual components of the binary system. The vi- 
sual extinction Ay and the efifective stellar temperature used 
in the model are within the 1 cr-deviations (1.0 mag and 400 K) 
of the previously derived results of Prato et al. (2003). The 
NIR flux that results from our model deviates by 60% from the 
measurements (s. Fig. |6]l. The accretion luminosity is equal to 

^acc ~ Lq. 
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Fig. 5. SED and MIR visibility for the projected baselines of B = 59 m, Z? = 65 m, and B = 74 m obtained from the observations 
and the model of HBC 639. We assume that the photometric measurements in the FIR range are ascribed to the infrared companion. 
The light gray curves/bars represent a model with a disk parameter/? =1.1 instead of using j8 = 1.0 (grey colored curves/bars). 



A second model for S CrA N (grey curve/bars in Fig.|6ll with 
j6 - 0.90 and /iioo = 12 AU does not reproduce the visibility 
measurements as well as the first model (deviations from the 
measurements are 6%) but the SED for A < 14 fj.m can be simu- 
lated more accurately. In contrast to HBC 639, the FIR flux must 
then originate in the Southern component. Furthermore, to re- 
produce the visual wavelength range, the second model still al- 
lows an inclination angle of up to ^ 45°, while only lower 
inclinations angles of§< 10° can be used in the first model. 



4.7. SCrAS 



The Southern component of S CrA is the fainter component in 
the NIR and MIR wavelength range. An investigation of the Bry 
line showed that this component is as active as the Northern 
component (Prato & Simon |1997| l. A value for the accretion 
rate had not yet been derived. The accretion luminosity in our 
model is Lace - 0.2 Lq. The corresponding accretion rate of 
M = 4 X lO^^Moyr"' is akeady an upper limit because higher 
values would produce an increase in the disk iiTadiation and a 
decrease in the visibilities. If accretion is not considered in our 
model, the MIR visibilities increase by only 2%. Figure |7]repre- 
sents our best-fit model. 



5. Radial gradient of the dust composition in 
circumstellar disks around TTauri stars 

For each single interferometric observation with MIDI and for 
each single telescope, one uncorrected N band spectrum can 
be obtained. For each interferometric observation, we also ob- 
tained a correlated spectrum reflecting the flux emitted by re- 
gions that are not spatially resolved by the interferometer. An in- 
creasing eff'ective baseline length B of the interferometer results 
in a higher resolution. The uncorrected, i.e., single-dish spectra 
as well as the coiTelated spectra contain spectral contributions 
of the N band from the entire disk, while the contributions from 
the hotter and brighter, i.e., inner regions increase with increas- 
ing effective baseline length B assuming a homogeneous, axial- 
symmetric disk. It was shown by Schegerer et al. (2008 ) that the 
local silicate dust composition in the circumstellar disk around 
RY Tau depends on the radial distance from the star, i.e., the rel- 
ative contribution of crystalline and large dust grains increases 
towards the star|f| Small and large dust grains are assumed to 
have a radius of 0.1 yum and 1.5 fim, respectively. In the follow- 
ing, we investige whether a corresponding correlation between 



' In this study, crystalline and large grains are called evolved dust 
grains. 
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Fig. 6. SED and MIR visibilities for the projected baselines of B = 56 m, and B = 70 m obtained from the observations and the 
model of S CrA N (dark grey line/colored bars). A second model with j3 = 0.90, and hioo = 12 AU (light gray lines/colored bars) 
can even improve the reproduction of the NIR fluxes. 
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Fig. 7. SED and MIR visibilities for a projected baseline of B = 72 m obtained from the measurements and our best model for 
SCrAS. 



radial location and dust composition can also be found in the 
TTauri stars of our sample. 

The (averaged) uncorrected spectra and the correlated spec- 
tra are shown in Fig. [8] All uncorrected and correlated spec- 
tra exhibit the silicate emission feature apart from the correlated 
spectra of HD 72106 and HBC 639 that were obtained with the 



longest interferometric baseline. The uncorrected spectra ob- 
tained during different nights differ in absolute flux scale by 
4% for SCrAN and up to 35% for GWOri, while the shape 
of the spectra does not change. We assume that the variation in 
the absolute flux scale of a spectrum is caused by an erroneous 
photometiic calibration during data reduction. Przygodda (i2004] l 
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indeed confirmed that the scale factor of the photometric cal- 
ibration can vary by more than 15% during a single night. 
If available, spectra acquired by the Thermal Infrared Multi 
Mode Insti-ument2 (TIMMI2) at ESO's observatory LaSilla 
(Przygodda et al. 120031 Schegerer et al. 2006 ) are also plotted in 
Figs.[8]and|9] The spectra obtained with TIMMI 2 can be used in 
the following comparison. With respect to the TIMMI 2 data, the 
shapes of the uncorrelated MIDI spectra are largely preserved 
for GWOri (deviation: 6%), RULup (9%), SCrAN (18%), and 
SCrAS (25%) and differ only in absolute scale. However, the 
shapes of the spectra of DR Tau, and of HBC 639 in particular, 
differ significantly in a way that could be due to an intrinsically 
temporal variation. We note that the TIMMI 2 spectra were ob- 
tained ~3 years before our MIDI observations. However, the dif- 
ference could also be caused by a doubling of the spatial res- 
olution power using a 8.2 m-single-dish of MIDI in contrast to 
the 3.6 m-telescope of TIMMI 2. Therefore, observations with a 
single-dish of the VLTI allows us to observe exclusively more 
central regions of the source, where more evolved grains are as- 
sumed to exist. We should mention that a spatial resolution of 
~43 AU at a distance of ~140pc can theoretically be achieved 
with a single-dish observation with the VLT. With respect to 
the TIMMI 2 spectrum, the flattening of the MIDI spectrum of 
HBC 639 could originate from an increase in the mass contribu- 
tion of large silicate grains in the more central regions consider- 
ing their specific absorption efficiency k (Dorschner et al. 1995 ). 
Our latter assumption is also confirmed by the correlated spec- 
tra of the objects HD72106B and HBC 639 obtained with the 
longest interferometric baseline. These spectra are flat and do 
not show any silicate emission band. We refer to a study of Min 
et al. ( I2OO6I 1, where it was found that the silicate band disappears 
if the averaged dust grain size exceeds ~4 jum. 

The origin of the evolved dust contributing to the TIMMI 2 
spectra can be derived by considering the following compari- 
son with the MIDI spectra. A black body B{v, T) that represents 
the putative underlying continuum (Schegerer et al. 120061) and 
is fitted to the extended wings of the emission features F(v) of 
TIMMI 2 as well as correlated and uncorrelated MIDI spectra, is 
subtracted from N band. We derive normalized features using: 



„(v) = 1 + 



(F(v) - B(v, T)) 
B(v, T) • 



(4) 



This normalization procedure preserves the shape of the emis- 
sion feature (Przygodda et al. I2003I I. According to Schegerer et 
al. (2006), the contributions of amorphous and crystalline dust 
to the TIMMI 2 spectra are determined using a fitting rou- 
tine, where the absorption efficiencies k of dust are linearly com- 
bined. In the context of the latter findings, we additionally sub- 
tract the contribution of small amorphous dust from the normal- 
ized TIMMI 2 spectra and compare the resulting spectra with 
the normalized uncorrelated and correlated MIDI spectra0 In 
Figs.[TO]and[TT] we only compare the normalized TIMMI 2 spec- 
tra with both the normalized uncorrelated and normalized corre- 
lated MIDI spectra obtained with the longest effective baseline 
length. With respect to our results, all sources can be divided into 
three classes: 



originate therefore in more central regions. Since the normal- 
ized TIMMI 2 spectra similarly correspond to the normalized 
uncorrelated and normalized correlated MIDI spectra, a fur- 
ther evolution of dust in regions that MIDI is able to study, 
cannot be determined. DR Tau and S CrA N belong to this 
class. 

ii. The normalized TIMMI 2 and MIDI spectra of HBC 639 dif- 
fer at ~10.8//m, where the crystalline compound enstatite 
generally shows a specific emission feature. If the remain- 
ing feature at ~10.8 yumin the normalized TIMMI 2 spectrum 
can be ascribed to enstatite, our comparison suggests that en- 
statite has its origin in more outer disk regions. These outer 
regions can already be resolved with a single-dish of MIDI 
and do not appear in the MIDI spectra. S CrA S could also be 
ascribed to this class by considering its correlated, normal- 
ized spectra. This result is based on the assumption that the 
spectral resolution power of MIDI is high enough to resolve 
the ~10.8 |/m feature of enstatite only in the inner regions. 

iii. An increasing flattening of the emission feature and an in- 
creasing amplitude of the forsterite feature at 11.3 yum with 
increasing baseline length B suggest grain evolution towards 
more central regions (e.g., GWOri). The comparison with 
the coiTesponding normalized TIMMI 2 spectra shows, how- 
ever, that the contribution of small amorphous dust grains, 
emitting a triangular emission feature with the maximum at 
9.8 fim, is still large even in the correlated spectrum that were 
obtained with the longest baselines. GWOri, RULup, and 
S Cr A S belong to this class. 

We note that all of these findings are based on the assumption 
that the underlying continuum of the emission feature can be 
represented by a single black body. This procedure is circum- 
stantially discussed in Schegerer et al. ( I2006I I. 

Considering a measurement error in the visibility of cr w 0.1 
as well as the low signal-to-noise ratio of ~4 and the low spec- 
tral resolution reached, we avoid fitting the uncorrelated and 
coiTelated MIDI spectra using a linear combination of different- 
sized, amorphous and crystalline silicate grains. The latter pro- 
cedure was presented by Schegerer et al. (2006) and Schegerer 
et al. (2008). However, we conclude that the local silicate dust 
composition in the circumstellar disks around the TTauri ob- 
jects of our sample depends on the radial disk location. Features 
of crystalline and large silicate grains indeed originate predom- 
inantly in more central disk regions of several AUs in size. 
However, the features of enstatite could originate in more outer 
regions, and the spectral contribution of small, non-evolved dust 
grains from the inner disk region could remain high. 



6. Summary 

We have presented MIR visibilities of seven pre-main-sequence 
stars observed with MIDI at the VLTI. We modeled the SEDs 
and the spectrally resolved MIR visibilities of these YSOs, and 
in particular T Tauri objects. The density distributions of the cir- 
cumstellar disks were derived using the parametrized approach 
of Shakura & Sunyaev (119731 Eq. [TJ. The results of this study 



i. The normalized TIMMI 2 spectra reproduce well the normal- 
ized uncorrelated and normalized correlated MIDI spectra. 
The evolved grains that contribute to the TIMMI 2 spectra 



* For simplicity, the modified TIMMI 2 spectra after normalization 
and after subtraction of the contribution of small amorphous grains are 
called normalized TIMMI 2 spectra. 



i. All objects could be spatially resolved with MIDI. 

ii. We showed that the results of highly spatially resolved ob- 
servations in the MIR range probing small-scale structures 
in the range of several AUs, and the SED from the visual to 
the millimeter wavelength range can simultaneously be re- 
produced by a single model except for one source (see vi.). 
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Fig. 8. Uncorrelated and correlated N band spectra obtained from the objects that were observed with MIDI for this study (black 
lines). Correlated spectra are designated by the additional information about the baseline length B (Table|2]i. Because of the remain- 
ing ozone band in several spectra the wavelength interval from ~9.3 //m up to ~9.9 //m is cut and underlaid with gray color. The 
gray lines represent data formerly obtained with TIMMI2 (Przygodda et al. 120031 Schegerer et al. l2006l l. 



We note that the SED does not generally provide any spa- 
tial information. The number of models that simulate the 
SED, solely, can be strongly reduced by the measurements 
with MIDI. Furthermore, the results of previous investiga- 
tions based on large-scale observations in the range of sev- 
eral 10 AUs, were confirmed by our modeling results in most 
cases. Any differences between the results of previous mea- 
surements and our modeling approach (e.g., the accretion 



rate of RU Lup) could be caused by a time-dependent varia- 
bihty. 

iii. For five of seven objects, the modeling approach of a purely 
passive disk is insufficient to reproduce the SED and the MIR 
visibilities of the sources. The implementation of accretion 
effects also significantly improves the simulation of the mea- 
surements. 
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Fig. 9. Continuation of Fig. [8] 



iv. For two objects in this study, i.e., HD 72106 B and HBC 639, 
the approach of a passive disk without accretion is sufficient. 
HD 72106 B is 10 million years old and its low accretion rate 
and high mass contributions from crystalline and large dust 
grains, have already been mentioned in previous investiga- 
tions (Vieu-a et al. 120031 Schutz et al. 120051 ). In contrast, vii. 
HBC 639 is still a relatively young object (Gras-Velazquez 
& Ray al. l2005l l but belongs to the class of weak-line T Tauri 
stars. 

V. For some targets (e.g., DRTau, and S CrAN), the MIR data viii. 
do not constrain all model parameters. However, further 
modeling constraints can be obtained from additional mea- 
surements, e.g., NIR visibilities as shown for our model- 
ing approach for DR Tau, where PTI data were additionally 
taken into account. 

vi. The modeling result for the source GWOri illustrates the 
need for an individual extension to our approach. According 
to Mathieu et al. (1997), GWOri has a stellar companion 
at a (projected) distance of ~1 AU. Therefore, in a second ix. 
model we considered a combination of a circumstellar and 
circumbinary disk. The circumstellar and circumbinary disk 
are determined by the same disk parameters j3 and /iioo and 



are seperated by a dust-free gap. However, a model that re- 
produces all the measurements could not be found. We rec- 
ommend acquiring additional high spatial resolution obser- 
vations including phase measurements to clarify the complex 
inner-disk structure of this object. 

Our measurements of RU Lup can be simulated equally with 
an active disk model with or without a dust-free gap. A 
disk model with a gap was formerly proposed by Takami et 
al. ( 120011 ) based on spectro-astrometric observations. 
Each single component of the binary system S CrA could 
be seperately observed with MIDI. Considering inner disk 
regions on small scales that could be resolved with MIDI, 
basic differences between this binary system and other ob- 
jects in our sample without a companion could not be found. 
Only the SEDs of the sources HD72106B, HBC 639, and 
S CrAN that decline strongly at A > lOyum support the idea 
that a stellar companion truncates the outer disk regions of 
this sources. 

The relative mass contribution of evolved dust in the systems 
of fainter T Tauri objects is higher in the inner disk regions 
close to the central star than in the outer regions. However, 
enstatite could already be enriched and stimulated in outer 
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Fig. 10. Uncorrelated and correlated MIDI spectra after subtraction of a black body that represents the underlying continuum and is 
fitted to the extended wings of the emission features (black lines). The gray curves represent TIMMI 2 data (Przygodda et al. 120031 ) 
after subtraction of a black body and the spectral contribution of small amorphous dust grains found in the preceeding study of 
Schegerer et al. (|2006|). For clarity, we only show the results obtained for the single-dish spectra and the correlated spectra that were 
obtained for the longest effective baseline length. The vertical lines designate features of enstatite (at 9.2 fj.m, 10.8 jum, and 11.1 yum) 
and forsterite (at 1 1.3 fim). 



disk regions and the spectral contribution from non-evolved 
dust grains could still be high in the inner regions. 
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Appendix A: Previous measurements: 

DRTau: The 2.5 million-year-old object DRTau belongs to the Taurus-Auriga star-forming region at a distance of ~140pc (Siess 
et al. 1 19991 1 ■ Numerous publications illustrate that DRTau is one of the most well-studied TTauri objects. Strong veiling has been 
found in the visual and NIR wavelength range: the flux in V band, for instance, exceeds the intrinsic stellar flux by a factor of five 
(Edwards et al. |19931 l. Apart from the excessive flux in the visual, the profiles of several emission lines, such as the Pfy-line, points 
to material accreting onto the central star while the profiles of further emission lines such as the He I line (Kwan et al. ,2007) and 
the Ha lines (Vink et al. I2005I I are evidence of outflowing stellar/disk winds. Based on a modeling study, Edwards et al. (I1993I I 
determined an accretion rate of M = 7.9 x lO^^Moyr"', while a mass of Mwind = 2.5 x lO^'^Moyr"' is lost by stellar/disk winds. 
A strongly collimated outflow was also found by Kwan & Tademaru ( 198 8). DRTau is photometrically and spectrally variable on 
short terms, i.e., in the range of weeks (AV =1.3 mag; Grankin et al. 120071 Eiroa et al. 20021 Smith et al. |1999l l. This variability is 
also ascribed to the formation and movement of stellar spots on this star (Ultchin et al.^l997|l. 

GWOri: GWOri, also known as HD 244138, belongs to the star-forming region B 30 in a ring-shaped molecular cloud close to 
A Ori. According to Dolan & Mathieu (1200 11 1, the ring shape has its origin in a central supernova explosion 1 million years, ago. As 
GW Ori is ~1 million years old (Mathieu et al. |199H ), the supernova explosion could cause the formation of this object. Considering 
a stellar luminosity of several tens of solar luminosities (Mathieu et al. 119971 Calvet et al. I2004I I. GW Ori is one of the most 
luminous YSOs with a spectral type of GO. Using theoretical evolutionary tracks, a stellar mass of 2.5 M© (Mathieu et al. |1991| l 
and 3.7 M© (Calvet et al. 2004) could be derived. GWOri is a spectroscopic binary (Mathieu et al. 1991 ). The companion with a 
mass of 0.5 M© up to 1 M© orbits the primary at a projected distance of ~1 AU in 242 days. Mathieu et al. ( 1 19971 ) used two different 
modeling approaches to reproduce the SED of this system. In their first approach, the secondary creates a (gas and dust free) gap 
between 0.17 AU and 3.3 AU. Their second modeling approach, where the circumbinary disk was replaced by a spherical envelope, 
disagreed with subsequent millimeter measurements at the James Clerk Maxwell Telescope (Mathieu et al. 1995 ). A disk mass of 
0.3 Mq could be derived using the latter set of millimeter measurements, where the object could be spatially resolved. An outer 
radius of 500 AU and an inclination angle of§~n° were determined. Artymowicz & Lubow ( 1 19941 ) showed that the disk gap of 
GW Ori proposed by Mathieu et al. (119911 ) cannot be explained by a theoretical modeling study of tidal forces. Calvet et al. (120041 ) 
measured an accretion rate of M > 2.5 x 10"^ Moa"', but the source does not reveal any veiling in the visual. To stabilize a high 
accretion rate for several 100, 000 years, Gullbring et al. (2000) proposed the existence of an additional massive envelope where an 
inner cavity enables the observation of the inner disk edge. GWOri is only weakly variable (AV - 0.2 mag; Grankin et al. 120071 ). 
Mathieu et al. (119911 ) pointed to a second companion with a period of 1 , 000 days that was found by the movement of the center of 
gravity in the system. 

HD72106B: The obje ct HD 72106 is a vis ual bin ary (angular distance 0.78" , PA = 199.8°) in the Gum nebula at a distance of 
288!^|J^pc (Torres et al. [T995l Hartkopf et al. [T996l Fabricius & Makai-ov|2000l). The Ha -emission as well as the infrared excess are 
ascribed to the visually fainter B-component (AV = 0.8 mag), while the A component already belongs to the main sequence (Vieira 
et al. 2003 ; Wade et al. 2005 ) and shows a strong magnetic field that was found with spectropolarimetry (Wade et al. 2005 ; Wade et 
al.|2007 ). The faint Ha emission line as well as the broad 10//m silicate emission band that can be compared with the silicate band 
of the comets Hale-Bopp and Halley point to the advanced evolutionary status of the B component as a YSO (Vieira et al. 120031 
Schiitz et al. 20051. In particular, Schiitz et al. (120051 ) found larger amounts of enstatite (~50%) that eff'ectively contributes to the 
silicate band. This large spectral contribution of enstatite has only been found for the evolved Herbig Ae/Be objects, HD 100546 
and HD 179218, so far. Folsom (200 7 ) intensively studied this binary system. 

RULup: RULup is a classical TTauri star in the star-forming region Lupus. Visual and millimeter measurements showed that the 
object does not have a (remaining) circumstellar envelope (Giovannelli et al. '1995 ; Lommen et al. 120071 ). A secondary could not 
be found with speckle interferometry down to a minimal distance of 0.1" in the NIR range (Ghez et al. 1 19971 ) and with the Hubble 
Space Telescope (Bernacca et al. 119981 ). Broad visual emission lines as well as flux variations in the U and V band are evidence of 
accreting material. Lamzin et al. (119961 ) determined a mass accretion rate of M = 3 x 10"^ Moyr The stellar magnetic field of 
the object has a strength of 3 kG (Stempels et al. 120071 ). Absorption lines shifted to longer wavelengths results from an outflowing 
stellar wind (Herczeg et al. 12005] ). 

HBC639: HBC639, also known as DoAr24E, belongs to the star-forming region p Ophiuchi. It is a class II object (McCabe et 
al.|2006 ) with an age of between 1 and 3 million years (Gras- Velazquez & Ray 2005 ). Because the Ha line reveals an equivalent 
width of ~5 A, HBC 639 is a weak-line T Tauri object, where the circumstellar material has evolved more rapidly than in classical 
TTauri stars of similar age. However, the infrared excess points to a remaining circumstellar disk in the system. The accretion 
rate inferred from the width of the Pa/S and Bry lines is low: M = 6 x 10 '^Moyr ' (Natta et al. 120061 ). HBC 639 has an infrared 
companion at an angular distance of 2" (320 AU for 160 pc) and a position angle of PA = 150° measured from North to East 
(Reipurth & Zinnecker 1993). The brightness of the secondary increases in the infrared and exceeds the brightness of the primary 
already in L band (Prato et al. 120031) . Chelli et al. ( 119881 ) assumed that the primary effectively contributes only to the NIR and MIR 
flux of the system. Polarimetric measurements in the K band showed that both components have a circumstellar disk with almost 
identical position angles (PA ^ 12,5°; Jensen et al. 120041 ). The companion is a class I object and active (Prato et al. 120031 ). Using 
the speckle-interferometric technique in K band, Koresko (120021 ) found a second companion close to the secondary. Both latter 
companions have a similar brightness in K band. 

S CrA: The source S CrA belongs to the Southern Region of the Corona Australis Complex (e.g., Chini et al. 120031 ). The source was 
already defined as a TTauri object by Joy et al. (1945) who pointed to an infrared companion at a projected distance of 169 AU and 
at a position angle of PA - 149° (Joy & Biesbrock 1944; Reipurth & Zinnecker |1993l ). Highly spatially resolved observations in 
the NIR wavelength range showed that both objects have an active circumstellar disk (Prato et al. 120031) . The spectral fines of both 
objects indicate that they have similar shapes and depths. Both components are probably coeval (Takami et al. 120031 ). S CrA is a 
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Table A.l. Photometric flux measurements of DR Tau. 



wavelength [fim] 


flux [Jy] 


Ref. 


0.36 


0.038 ± 0.021 


1 


0.44 


0.050 ± 0.049 


1 


0.55 


0.083 ± 0.060 


1 


0.64 


0.13 ±0.06 


1 


0.79 


0.21 ± 0.01 


1 


1.25 


0.45 ± 0.01 


2 


1.65 


0.77 ± 0.04 


2 


2.20 


1.13 ±0.02 


2 


3.50 


2.07 ± 0.82 


1 


3.60 


1.86 ± 0.20 


7 


4.50 


1.89 ±0.15 


7 


4.80 


2.41 ± 1.10 


1 


5.80 


1.27 ± 0.01 


3 


8.00 


1.77 ±0.20 


7 


12.0 


3.16 ±0.25 


4 


25 


4.30 ± 0.29 


4 


60 


5.51 ±0.54 


4 


100 


6.98 ± 0.63 


4 


200 


3.98 ± 0.84 


5 


450 


2.38 ±0.17 


6 


729 


0.40 ± 0.08 


8 


850 


0.53 ±0.01 


6 


1300 


0.16 ±0.01 


6 



References - 1: Kenyon & Hartmann tl995i ; 2: 2 MASS catalogue (Cutri et al. |20031 ; 3: Haitmann et al. |2005h 4: Gezaii catalogue fT999t ; 5: ISO data archive; 6: Andrews & 

Williams (2005); 7: Robitaille et al. (2007); 8: Beckwith et al. fT990) 

Table A.2. Photometric flux measurements of GW Ori. The value that is marked with the symbol "J," is an upper flux limit. 



wavelength [fim] flux [Jy] Re: 



0.36 


0.068 ± 0.014 




0.44 


0.21 ± 0.02 




0.55 


0.43 ± 0.02 




0.70 


0.61 ± 0.03 




0.90 


0.82 ± 0.05 




1.25 


1.31 ±0.04 


2 


1.65 


1.47 ± 0.04 


2 


2.20 


1.47 ±0.04 


2 


3.50 


2.45 ± 0.09 


3 


4.80 


2.7 ± 1.2 


4 


18 


17 ±3 


5 


25 


20.5 ± 1.2 


6 


60 


31.5±4.1 


6 


100 


35i 


6 


350 


5.0 ±0.6 


7 


450 


3.5 ± 0.4 


7 


800 


0.9 ±0.1 


7 


850 


1.0 ±0.1 


7 


1100 


0.29 ± 0.03 


7 


1360 


0.26 ± 0.06 


7 



References- 1: Calvet et al. t2002) ; 2: 2MASS catalogue (Cutri et al. l2003l ; 3: Rvdgren & Vrba fT983l ; 4: Cohen & Schwartz (1975); 5: Cohen fT973) ; 6: IRAS catalogue IT985) ; 

7: Mathieu et al. fT995l 

YY Ori object, i.e., the emission lines are asymmetric and shifted to longer wavelengths. These lines arise from material that accretes 
onto the central star. The measured spectral variability is another hint of non-continuous accretion process. Prato & Simon (1 19971 1 
showed that only an infalling circumbinary envelope provides enough material for accretion in the long-term. 
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Table A.3. Photometric flux measurements of HD72106B. The value that is marked with the symbol "X" is an upper limit. 



wavelength \pm] 


flux [Jy] 


Ref. 


0.36 


0.80 


1 


0.44 


1.67 


1 


0.55 


1.45 


1 


0.79 


1.32 


1 


0.90 


0.89 


1 


1.23 


0.61 


3 


1.25 


0.67 + 0.01 


2 


1.63 


0.45 


3 


1.65 


0.50 ± 0.02 


2 


2.19 


0.41 


3 


2.20 


0.43 ± 0.01 


2 


3.78 


0.39 


3 


4.66 


0.44 


3 


8.36 


1.74 


3 


9.67 


1.45 


3 


12.0 


2.23 ±0.13 


3,4 


25 


3.62 ± 0.21 


3 


60 


1.88 ±0.11 


4 


100 


16.8^ 


4 



References - 1 : Torres et al. 119951 ^ 2: 2 MASS catalogue (Cutri et al. l2003l : 3: Gezari et al. (1999); 4: IRAS catalogue fT985l 

Table A.4. Photometric flux measurements of RU Lup. 



wavelength [/im] 


flux [Jy] 


Ref. 


0.36 


0.016 


1 


0.44 


0.073 


1 


0.55 


0.115 


1 


0.64 


0.184 


1 


0.79 


0.290 


1 


1.25 


0.51 ±0.01 


2 


1.25 


0.72 ± 0.02 


3 


1.25 


0.41 ±0.01 


3 


1.65 


0.76 ± 0.03 


2 


1.65 


1.02 ± 0.03 


3 


1.65 


0.62 ± 0.02 


3 


2.20 


0.89 ± 0.02 


2 


2.20 


1.29 ± 0.02 


3 


2.20 


0.75 ± 0.03 


3 


3.50 


1.16 


4 


3.50 


1.77 ±0.07 


3 


3.50 


1.0 ±0.10 


3 


4.80 


1.5 


4 


4.80 


1.5 ±0.2 


3 


4.80 


0.6 ± 0.3 


3 


15 


2.22 


5 


20 


2.80 


5 


25 


4.64 


4 


60 


4.68 


4 


100 


5.70 


4 


1300 


0.197 ±0.007 


6 


1400 


0.159 ±0.010 


7 



References- 1: Gahm et al. 119931 : 2: 2MASS catalogue tCutri et al. l2003< ; 3: Giovannelli et al. 119951 : 4: Gezari catalogue |1999» : 5: Gras- Velazquez & Ray (2005); 6; Nurnberger 

et al. fT997l ; 7: Lommen et al. (2007) 
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Table A.5. Photometric flux measurement of HBC 639. The flux can be ascribed to the main component up to wavelengths in L 
band. 



wavelength [/im] flux [Jy] Ref. 

036 ooooi r 

0.44 0.0008 1 

0.55 0.0053 1 

0.64 0.018 2 

0.79 0.065 2 

1.25 0.33 ±0.10 3 

1.65 0.77 ±0.16 3 

2.20 0.92 ± 0.09 3 

3.50 0.75 ± 0.07 4 

4.80 0.98 5 

7.70 0.60 5 

10.0 1.06 5 

12.0 1.87 5 

12.8 2.57 5 

15.0 2.41 5 

20.0 2.30 5 

25 5.91 6 

60 35 6 

100 40 6 

800 0.037 ± 0.008 7 

1300 0.065 ± 0.019 8 



References- 1 : Bouvier & Appenzeller 119921 : 2: Rvdgren 119761 : 3: Prato et al. 120031 ; 4: Mc Cabe et al. 120061 : 5: Gras- Velazquez & Ray 120051 : 6: Gezari Catalogue fT999l : 

7: Jensen et al. (1996); 8: Nurnbeiger et al. {19971 



Table A.6. Photometric flux measurement of S CrA N and S CrA S. For the photometric observations of wavelengths A < 0.64 yum 
and A> 12 /im, both components could not be spatially resolved. 



wavelength \jim] 


SCrAN SCrAS 
flux [Jy] flux [Jy] 


Ref. 


0.44 


0.0089 


1 


0.55 


0.13 


2 


0.64 


0.09 


1 


1.25 


0.57 ± 0.02 0.28 ± 0.02 


3 


1.65 


1.02 ± 0.04 0.50 ± 0.03 


3 


2.22 


2.44 ±0.21 0.74 ±0.09 


4 


3.50 


3.99 ±0.23 1.15 ±0.03 


4 


8.80 


2.90 ±0.20 1.59 ±0.17 


4 


10.6 


3.19 ±0.20 1.27 ±0.06 


4 


12.0 


3.97 


5 


12.0 


4.88 ±0.59 


5 


12.5 


1.50 ±0.10 


5 


25 


9.21 ±0.64 


5 


25 


8.52 


5 


60 


1.81 ±0.16 


5 


60 


2.7 ± 4.5 


6 


60 


3.32 ±3.1 


6 


90 


2.87 ± 2.85 


6 


100 


2.79 ± 0.75 


5 


170 


8.29 ± 22.24 


6 


450 


2.75 ±0.38 


7 


800 


0.77 ± 0.05 


7 


1100 


0.36 ± 0.03 


7 



References- 1: Denis data base (2005) ; 2: Takami et al. 120031 : 3: Prato et al. 120031 ; 4: McCabe et al. 120061 ; 5: IRAS catalogue 119851 : 6: ISO data archive; 7: Jensen et al. 119961 



